Universality Classes of Diagonal Quantum Spin Ladders 
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We find the classification of diagonal spin ladders depending on a characteristic integer Np in 
terms of ferrimagnetic, gapped and critical phases. We use the finite algorithm DMRG, non-linear 
sigma model and bosonization techniques to prove our results. We find stoichiometric contents in 
cuprate Cu02 planes that allow for the existence of weakly interacting diagonal ladders. 
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The determination of the ground state and low lying 
excited states properties in strongly correlated systems is 
known to be a central problem in condensed matter since 
they are responsible for their thermodynamical proper- 
ties at low temperature, their interaction with weak ex- 
ternal fields and generically, for the presence or not of 
long-range correlations in the system. Thus, any time 
a new quantum system is brought about, their low en- 
ergy properties are one of the first issues at stake. In 
this work we adress this problem for the case of diagonal 
quantum spin ladders with antiferromagnetic Heisenberg 
interactions: 



iv/ The necklace diagonal ladder is the simplest mem- 
ber of the class having Np = 1 plaquettes. It is realized 
in both organic and inorganic compounds (e.g.: metal- 
free poly(m-aniline) or poly[l,4-bis-(2,2j6,6-tetramethyl 
l-4-oxy-4-piperidyl-l-oxyl)-butadiyne]) p. v/ Standard 
square ladders have been sinthesized in cuprate super- 
conductors like Srn-iCun+i02n,n S {3, 5,7, ••■} [^. 
Here we shall provide an explicit experimental proposal 
to study the whole class of diagonal ladders in cuprates. 

In this paper we cast the main results of our investiga- 
tions in a diagram of universality classes for AF diagonal 
ladders depending on their number of plaquettes Np-. 
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FIG. 1: a) Diagonal ladder with Np — 4 plaquettes. b) The 
Np — 2 diagonal ladder as a decorated zig-zag ladder, c) 
Rectangular ladder with n; = 2 legs. 

an array of sites hosting quantum spins S = and 
links, characterized by their number of plaquettes Np as 
depicted in Fig. ^ On the contrary, standard square 
ladders are characterized by the number of legs n; 
The interest in this type of spin systems is manifold: 
i/ They represent an alternative route to reaching 2D 
quantum physics from ID systems ii/ They ap- 

pear related to the physics of stripes in cuprate materials 
{e.g.Lai,Q-xNdoASrxCuOi) 0. iii/ They also show up 
in numerical simulations of the t — J model in planes [j| . 
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Odd Np Ferrimag., Gapless (F), Gapped (AF) 
Np — 2 {mod 4) AF, Haldane phase 



Even N„ 



Np ^0 {mod 4) AF, Gapless 
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The first classification of diagonal ladders occurs de- 
pending on whether they show ferrimagnetism or not: for 
Np odd they show ferrimagnetic order while for Np even, 
they do not. The reason is as follows: although all di- 
agonal ladders are bipartite systems, the number of sites 
in each sublattice is different. This unbalance between 
the number of sites in each sublattice is responsible for 
the Np odd ladders exhibiting ferrimagnetism. Moreover, 
these Np odd ladders present two types of excitations: 
a/ ferromagnetic gapless excitations which are the Gold- 
stone bosons associated to the ferromagnetic order of the 
ground state; b/ antiferromagnetic gapped excitations. 

The second classification corresponds to Np even lad- 
ders, which are not ferrimagnetic. Their classification 
is not as simple as the Np odd case since the gapless 
and gapped excitations do not coexists. In fact, for 
Np = {mod 4) the diagonal ladders are gapless sys- 
tems with long-range interactions and critical, while for 
Np = 2 {mod 4) they are gapped, with finite correlation 
length and we can describe them as a Haldane phase. 
Both types of excitations are antiferromagnetic. 

Therefore, it is apparent that the classification of di- 
agonal ladders is richer than that of standard square lad- 
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FIG. 2: The gap A for the Np — 2 diagonal ladder vs. cou- 
pling constant. Inset: idem for the string order paramteter. 
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FIG. 3: The gap A for the Np = 4 diagonal ladder vs. the 
length 1/L (unit cell) for J' — J, and Np = 2 for comparison. 



ders, which show gapless behaviour for ti; odd legs and 
gapped for rt; even In fact, we want to stress that the 
study of diagonal ladders is much more difficult than for 
square ladders. To see this, we introduce a topological 
deformation of the diagonal ladders as shown in Fig.^^) 
for Np = 2. This is achieved by flattening the outer links. 
The net result is a new sort of zigzag ladder that we call 
decorated zigzag ladder due to the extra sites. In the 
horizontal legs we put coupling constants J while in the 
zigzag inter- leg links they are J'. This distribution must 
be compared with the corresponding distribution for the 
square ladder in Fig. ^) where J' stands for the rung 
coupling constants. The difficulties in studying diagonal 
ladders can now be attributed to the absence of a simple 
strong coupling limit J' / J :s> 1: when J :s> J' the ladder 
is equivalent to 2 weakly coupled chains (gapless point) 
similar to what happens with the square ladder; when 
J' ;» J it is again equivalent to one weakly coupled chain 
(gapless point) unlike the square ladder which has a well- 
defined strong coupling limit around singlet states on the 
vertical rungs producing a gapped phase. The absence of 
a strong coupling limit for diagonal ladders have impor- 
tant consequences for their study: neither it is possible to 
perform a perturbation theory as in 7| nor a mean-field 
calculation based on the condensation of spin singlets on 
the rungs as in j^] . All this makes diagonal ladders a hard 
problem compared to square ladders. One of the funda- 
mental open questions we address is: are there gapless 
phases around the limits J ^ J' and J' ^ J? We need 
hereby to resort to a series of numerical an analytical 
nonperturbative methods. 

DMRG: We focus on Np even diagonal ladders and start 
the DMRG calculations [13, EJ for the Np ^ 2 case 
which is equivalent to a decorated zigzag ladder Fig.^). 
We use the finite algorithm version of the DMRG method 
which is more accurate than the infinite version. We have 
devised a sweeping process through the lattice adapted 
to the topology of the diagonal ladders for arbitrary 



Np which improves the accuracy of the algorithm: (a) 
all sites should be included just once (i.e., a Hamilto- 
nian graph); (b) ideally, all links should be included 
also (i.e., an Eulerian graph), but since this latter is im- 
possible, we reduce our pretensions to (c) most of the 
links should be included in the path, favouring the ones 
in the "short direction" , (d) path links which are not 
graph links are allowed, but not favoured. The maxi- 
mum number of sites is = 122 (which corresponds 
to a horizontal length of L = 30 unit cells), the num- 
ber of states kept is m = 200, discarded weights of the 
density matrix are always below 10~® and we achieve a 
convergence of 6 digits in the energies of the eigenstates. 
Moreover, we have used an independent Lanczos method 
to check our DMRG results for small number of sites up 
to iV = 28. In Fig. H we present our DMRG results for 
the gap A = E{L, S ^ 1) - E{L, 5 = 0) between the 
first excitated state (spin triplet) and the ground state 
(spin singlet) as a function of the coupling constant ra- 
tio J/ J'. Each point in the plot is the result of a ex- 
trapolated fit of the gap as a function of 1/L given by 
A(L) = oo + aiL~^ + a2L~^. With these results we can 
answer now the above open question: we clearly see that 
the gap vanishes only at the extreme cases of J' ^ 
and J — > 0. Thus, we have a dense gapped phase VJ'/J, 
except at two limiting points. 

The next interesting fact we find is the existence of a 
maximum corresponding to the isotropic point J = J' . 
This must be compared with the results for the standard 
zigzag ladder in where the physics of the gapped 
phase is dominated by the Majumdar-Gosh point (fully 
dimerized state) and in addition, there is a gapless phase 
separated by a critical point (J/J')c ~ 0.24. The rea- 
son for the standard zigzag ladder having two phases is 
because it is a frustrated lattice, unlike our decorated 
zigzag ladder which is bipartite and thus nonfrustrated. 
This makes a big difference and is responsible for the 
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FIG. 4: Block-decomposition suitable for diagonal ladders. 



presence of an single gapped phase in our ladder. Like- 
wise, we may wonder whether our gapped phase in Fig. [21 
is a real Haldane phase |3 Q| or dimerized. To clarify 
this issue we have computed with DMRG the dimerized 
order parameter D"{L) for a given length as the differ- 
ence between two consecutive spin-spin correlators in the 
measured in the middle of the ladder as 



(3) 



where a — u,d,z stands for up, down or zigzag di- 
rections in the ladder. For the Majumdar-Gosh state 
this parameter takes a maximum value and detects its 
dimerization. In our Np = 2 ladder we find values 
D"''^{L = 25) w 4 • 10"^ D^(L = 25) « 2 • IQ-^. Thus, 
our ladders are not dimerized, a result compatible with 
a Haldane phase. Moreover, we have also computed with 
DMRG the string order parameter (SOP) defined as 



(5) 



where is an 0(3) sigma field with the constraint 1 P — 
1 , £ is the angular momentum, the spin wave velocity, 
g the counpling constant and 9 is the topological term 
which is defined only mod 2tt. These {va-,g,0) are phe- 
nomenological parameters that will depend on the mi- 
croscopic parameters like (J, J', S) via the sigma model 
mapping The behaviour of the system depends on 
whether the topological constant 9 = (mod 27r) for 
gapped or 6* = TT (mod 27r) for gapless. This mapping 
can be done with a real-space blocking procedure 
For the diagonal ladders we need to use the topologi- 
cal deformation introduced earlier in Fig.^) in order to 
bring them to the form of generalized decorated ladders 



with ni 



+ 1 legs as shown also in Fig. 0] where we 



also show the splitting of the original Hamiltonian 
into a block Hb and interblock Hbb parts. In order to 
make an appropriate ansatz we need to take into account 
several subtleties occurring in diagonal ladders: a/ the 
Neel state has AF character in the horizontal direction 
but Ferromagnetic in the vertical direction despite being 
all the couphngs J, J' positive (Fig.0J; b/ the even-site 
spin operators are inequivalent to the odd-site spins due 
to the peculiar zigzag interleg couplings (Fig. 0J. With 
these provisos in mind, we find the correct ansatz to be: 



S{L) := ((85* + S^)e"'^i=2'(^^'+^^''(S^ + S^)). (4) 

Non- vanishing values of l@J detect the hidden order char- 
acteristic of the Haldane phase. In Fig. 13 (inset) we plot 
the extrapolated values S{oo) and find finite values in 
whole range corresponding to the gapped phase. It is re- 
markable that our values of the SOP are bigger than for 
the square ladder, while the gap in the diagonal ladder is 
more than ten times smaller than in the square ladder. 

For Np = 4 diagonal ladder, we have performed similar 
DMRG calculations up to = 184 sites with m = 400 
states. Our results are shown in Fig. |2| where we present 
the corresponding gap vs. 1/L, and also the same for 
Np = 2 for comparison. We find that the extrapolated 
value points to a negligible gap for the Np = 4 ladder. 
This is a crucial difference w.r.t. the previous case and it 
means that there is a sharp difference between the whole 
class of Np even diagonal ladders. Is it possible to give 
an explanation for this difference? 

Non-Linear Sigma Model: We turn to the sigma model 
(NLSM) [i3 to answer this question. The low energy 
properties of a Heisenberg lattice can be mapped onto 
the effective Hamiltonian of a non-linear sigma model 



Sa{2n) = A^eix) + ea{x) + S {(j){x) + Mx)) 
Sa{2n + 1) = A°J{x) + laix) - S i(P{x) + Mx)) 

(6) 

for the spin operators within each block in terms of sigma 
model fields and extra massive fields (j)a , ta needed to ac- 
count for the correct number of degrees of freedom in 
each block. A^^° are amplidutes computed from spin 
wave theory in terms of the microscopic parameters. To 
fulfill the canonical commutation relations between the 
sigma model fields we need to impose the constraints 

Ea^a = Y.aK = 1, Y.a^a = 0, Ea = 0. After a 

lengthy calculation we find the result: 

B = 2^s(^^^. (7) 

This is independent of J and J', a rather remarkable 
result having in mind the differences introduced by the 
zigzag couplings. We also calculate the values of g and 
Va in Table The coupling constant has very interest- 
ing information that clarifies the issue of lack of strong 
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TABLE L NLSM data {va,g,9) in TVp-even diagonal ladders. 



coupling in diagonal ladders, namely: 



[const. 



Np even diagonal ladders, 

ni even square ladders, 
ni odd square ladders. 



(8) 

Thus, it vanishes in the strong coupling limit J' ^ J 
meaning that the system becomes less disordered. This 
is in sharp contrast with the behaviour of square ladders. 
With this NLSM calculations we see that Np even ladders 
can have a very different behaviour as shown in Fig.|21and 
the DMRG calculations in Fig. El 

Bosonization: We can obtain additional information 
about the gap formation in the Np = 2 diagonal lad- 
der (decorated zigzag) using abelian bosonization tech- 
niques. Each leg of the decorated zigzag ladder (Fig.^)) 
can be described in the continuum limit in terms of mass- 
less Bosonic fields Lpa{x),a — 1,2 with a non- interacting 
Hamiltonian Hq = ^ J2a=i,2 I [^li^) + {dxfaix))^] ■ 
The corresponding SU(2) currents 3i^'^\x) and stag- 
gered field na{x) satisfy Sa{x) = Jai^) + -^ai^) + 
(— l)"n(a;). For weak coupling regime J' /J ^ 1, we may 
switch on the interchain coupling 

Hint = J'^ [Si(2n+ 1) + Si(2n + 3)] -82(271-1-2), (9) 



and bosonize it with the result Hint ~ — 4 J'oo / dx ni (x) ■ 
n2(x), where we have kept only the strongly relevant 
staggered parts and the marginal current-current terms 
are discarded. The nice feature of this result is that when 
we compare it with the bosonization of the standard 
square 2-leg ladder Hint ~ J'o^o J dx ni(a;)-n2(x) 0| and 
the standard zig-zag ladder iJint « J'ao J dx [—{dni{x))- 
n2(a;) + (9n2(x))-ni(x)] we find that our decorated zigzag 
ladder {Np = 2) produce a gap of the same type as the 
square ladder (both unfrustrated) and not related to the 
standard zigzag ladder which is dimerized. Therefore, 
this supports the existence of a Haldane phase in our 
Np — 2 diagonal ladder. 

Experimental Proposal: Let us now turn to the experi- 
mental realization of diagonal ladders mentioned at the 
beginning based on cuprate materials. We find that it is 
possible to modify the basic CM02-stoichiometry of the 
cuprate planes in order to arrange the Cu sites and the 
O sites in links of 180 degrees with the precise shape of 



a diagonal ladder, while each diagonal ladder are weakly 
interacting among each other because they share links of 
Cu — O — Cu making 90 degrees that are much weaker. 
This construction for the Heisenberg interaction is 
based on the superexchange mechanism among Cu sites 
mediated by interstitial O sites which is stronger when 
the sites Cu—O are aligned in 180 degrees. An example of 
our construction for the first Np = 1,2 diagonal ladders 
^ is shown in Fig. corresponding to a stoichiometry 
of CU3O4 and CU2O3, respectively. 
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FIG. 5: Np — 1 (left) and Np — 2 (right) weakly interacting 
diagonal ladders in cuprate planes. Solid circles are Cu sites 
and the lattice sites are the O interstitial sites. 



The outcome of our general construction is summa- 
rized in the following series of cuprate materials based 
on Stroncium: 



Srn-lCUr, 
Srn-lCUr. 



-l02n 
-\Oln 



n-2,4,. 
n = 3, 5, ■ 



A^p = l,3, 
A^p = 2,4, 



(10) 

Interestingly enough, the series of cuprates with n odd 
has been studied experimentally as a set of square lad- 
ders with variable number of legs n/ , both even and 
odd. What we have found is that our diagonal ladders 
with Np even are allomorphic with all the square ladders. 
This makes sense since we have found both gapped and 
gapless behaviour in Np even diagonal ladders lO. It is 
remarkable that each of the series in our classification (j^J 
of universality classes corresponds to a different series in 
the same family of cuprate materials and it may allow for 
an experimental test of those complex behaviours found 
in diagonal quantum spin ladders. 
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